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TIME VARIABLE OSMOTIC PRESSURES PRODUCED BY 
COUPLED REACTIONS AS A POSSIBLE CAUSE ~ 
OF CELL DIVISION 


N. RASHEVSKY 


COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


It has been shown in a previous paper that coupled reactions may 
produce diffusion phenomena which are characterized by highly asym- 
metric distributions of concentrations in a spherical cell. Under certain 
conditions such diffusion fields are unstable, the asymmetries in concen- 
trations tending to increase indefinitely. This results in an increase of 
asymmetrically distributed osmotic pressures which may eventually re- 
sult in a division of a cell. Constriction without elongation, as in cleay- 
age, is studied. The process of division is brought about in this case by 
ordinary osmotic pressure, and not by the diffusion drag forces. 


In a recent paper (Rashevsky, 1948a) we have shown that two 
simple coupled reactions may produce time-variable diffusion fields 
in a spherical cell which may exhibit a highly asymmetric spatial 
structure. However, those fields are in general unstable and cannot 
last permanently. The only stable permanent distribution of concen- 
trations is a spherically symmetric one. 

Even temporary asymmetries of concentrations, when large 
enough, will produce an asymmetric distribution of ordinary osmotic 
pressures, and this may result in forces which will tend to deform 
the cell and perhaps even divide it into two. 

The purpose of the present paper is to discuss this possibility. 
The reader’s familiarity with the previous paper is presupposed. 

By adding instead of subtracting equations (8) and (9) of loc. 
cit., we find that 

WY; = Ci + Cie (1) 


satisfies the same differential equation as ®;. Hence, W; will be given 
by similar expressions as #;. If, contrary to the usual convention, 
we express ¢;, and ¢i2 in mol cm-* rather than in gm cmr, which does 
not change anything in the equations, the total osmotic pressure at 
each point will be proportional to ¥;. Hence the osmotic pressure p 
will be given as a function of time and coordinates by expressions 


similar to those which give 9; . 
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Separating p into a time-dependent part p* and a time-indepen- 


dent part 2 , pa oe eae 


p=P.+ Pp, (2) 
using the same notations as in loc. cit., we may write [loc. cit., equa 
tion (28) ]: ees: as al 

eras A (3) 


where »y is given by [loc. cit., equation (45) ]: 


gee [A? (Din oe Diz) i oe | y rE MDiDis ar 7 (4D is 
(4) 
i+ Ao2D ix) + Gy1Ae2'— Ay2%0, = 0. 
The quantity 4 is the eigenvalue, determined by equations (63) and 
(71) of loc. cit. For the particular case where Da = De = 0 and 
De = 0 (corresponding to 6, = 0 in loc. cit.), 4 is determined by 


1 


d ii k 
ik 4 _ 7) he yieaae ee hAT oe te (5) 
adr Vr m4 r=1p / T mh 


Without loss of generality we shall confine ourselves in this paper to 
that case. 
Consider the axially symmetric particular solution 


a 
p* = — Js,2(Ar) Po (cos 6)ev', (6) 
Ae 


which is a particular case of equation (52) of loc. cit. The osmotic 
pressure has two maxima near the poles (Figure 1, shaded area) and 


FIGURE 1 


a minimum along the equator HE’ (Figure 1). If for a given value 
of 4, say the smallest value, 2, , the quantity v is real and positive, 
then the difference in osmotic pressures between the polar and equa- 
torial regions will increase with time. By applying Betti’s equation 
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it can readily be shown that such a distribution of pressure as shown 
on Figure 1 will tend to elongate the cell along the line PP’ and con- 
strict it along the line EE’. 
The first question is whether positive values of » are possible. 
The affirmative answer to this is obtained by a simple example. 
Denote by x, the smallest positive root of Js(x) = 0 and by 
2,’ the smallest positive root of J'5.(“) = 0. When k, is very large 
the second term of (5) prevails, and J, is given by 
x 
WT fOr = (7) 
Yo 
When £ is very small or zero, the first term of (5) prevails and A, is 
given by the equation 


A i 

ae eS (Aro) ae seeae 

NOP PAV gee 
Equation (7) gives, with 7, ~ 10° cm, 4, ~ 103. Let us assume for 
a moment that the value of j,, as determined by (8) is also of 
the order of magnitude of 10°. The coefficients of both terms of (8) 
are then of the same order of magnitude (10+). The smallest roots 
of both equations 


J 5/2(ATo) =0. (8) 


J'so(x) =0 and Jss(x) =0 (9) 


are of the order of magnitude of unity. Since both J’s,. and Js. are 
oscillating functions, the smallest value of Ar, for which (8) holds 
will also be of the order of unity. 

Thus both for small and large values of k,, the value 4, is of 
the order of 10°. 

Now let us take as plausible values (Rashevsky, 1948b) Di; = Diz 
~ 10-7 cm? sec, 1, = Ais = Ae, ~ 107 sec, Geo = 0. With these val- 
ues and with 4, ~ 10%, we find for the positive root of (4) the value 
y ~ 107 sec. With such a value of » the pressure difference between 
polar and equatorial regions will increase about 3 times every 10 sec- 
onds. Smaller or larger values of » may still be obtained within the 
plausible range of values of the other constants. 

The next question is how to calculate the effects of such a vari- 
able pressure on the cell. As the cell deforms under the action of this 
pressure, the solutions which were derived for a spherical cell will no 
longer hold. We must resort to the approximation method. Its gen- 
eralization to problems of this type has been outlined recently (Ra- 
shevsky, 1948c). We may consider an oblong cell, divided into “‘polar” 
and “equatorial” regions as shown on Figure 2. To such a cell we may 
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apply the approximation method in the form given in the paper men- 
tioned above (Rashevsky, 1948c). We thus will obtain the variation 
of concentrations and hence of pressures in both regions as functions 
of the length and width of the cell. Then, assuming in the usual way 
that the deformation is slow, we can apply Betti’s equation. 


FIGURE 2 


Leaving such a study for a later time we shall here consider, as 
an illustration only, the case in which the cell is prevented from elon- 
gating either by a sufficiently rigid outer membrane or by other cells 
in a tissue. Even in the absence of elongation there will be a tendency 
toward constriction in the equatorial region HH’, if the osmotic pres- 
sure here inside of the cell becomes sufficiently low. Let us derive the 
equation of the constriction assuming, as a very crude approximation, 
that this constriction does not too appreciably affect the distribution 
of the concentration. 

Such an assumption may at least partly be justified by the con- 
sideration that along the plane HEH’, op/oé = 0. The gradient of the 
concentration is therefore zero in the direction of the normal to EE’. 
Hence an introduction of an impermeable partition along HE’ should 
not affect the diffusion field. A narrow constriction such as shown on 
Figure 3 may be considered as roughly equivalent to such a partition. 


Vas 


FIGURE 3 


(The mechanical analogue to this is the fastening of a vibrating sys- 
tem at a node or along a nodal line or plane.) 
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Consider the “neck” (Figure 3) as a cylinder of length 24 and 
a radius r. Such a distribution of pressures inside of. the cell as dis- 
cussed above will result in a net pressure — p, (directed inward) as 
applied to the “ends” of the “neck”, and a pressure — p, applied to 
the “sides”, with p, and p. both positive and p, < p.. Now applying 
Betti’s equation to the neck, we find (Rashevsky, 1948b, p. 142), 


1 dA 1 


——— = —— |, 10 
Adt  3nV : oy 
where V is the volume of the neck, 
V=2nrA, (11) 


and I, is given byt 
I,=V (p2—Di1) =VPD; P=Pe—Di> 0. (12) 


To the osmotic pressure we also must add the capillary pressure 
acting on the “sides”. That pressure is composed of a negative. part, 
— y/r , and of a positive part y/A. The latter is due to the fact that 
the curvature in the meridional plane of the neck is negative and its 
radius is of the order of magnitude of 4. Therefore the total capil- 


1 ae 
lary pressure on the “sides” is — »( = --) . This adds to J, a term 


Uf 
(t-3), w 


If the “neck” elongates at constant volume, then (Rashevsky, 
1948b, p. 149) 


-—=———, (14) 
rdt 2A dt 


and therefore altogether we have, combining (10), (12), (18), and 
14}: 
dr i v 
eee (a |e (15) 
dt 


Actually in a pure constriction such as considered here, the neck 
does not elongate, and the volume is not constant. Since the elonga- 
tion occurs because of a net lateral pressure, we may consider the case 


+Notice the difference in the expression for the volume here (277r,7,2) and 


4 
in N. Rashevsky, 1948b, p. 142 C WIT a2) 
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is equivalent to one in which a plastic rod elongates under lateral 
pressure, but as it elongates it is gradually chopped off at the ends 
so as to keep the length constant. Approximately equation (15) will 
still hoid. 

Equation (15) is readily integrated if we consider the case in 
which p varies very slowly so that it does not change appreciably dur- 
ing the process of constriction. We may then consider p as constant. 


Putting 
: ( "=a; 7 (16) 
67 7. A 3 67 : 
and integrating (15) with the initial condition 
P= 7p tore b=0s (17) 
we find 
aro + b b 
———— ee (18) 
a a 


Constriction occurs only if a > 0, or p — y/A > 0. The division 
ends when r = 0, which happens at 


a eI 
a og i 0. (19) 


In this case the constriction follows an entirely different course 
than when the cell is allowed to elongate freely (Rashevsky, 1948b, 
p. 148 ff.). Data on this type of constriction are not available but 
could readily be obtained in a similar manner as those for freely elon- 
gating cells. 

If p « e’' we obtain more complicated expressions. 

If » is real and positive only for 4 = 4,, and either negative or 
complex with a negative real part for all other 4’s, then the numer- 
ous other possible configurations which correspond to the individual 
terms of the general expression (65) of loc. cit. will exist for a while 
but gradually will disappear. As seen from (4), the damping is foe 
tive for sufficiently large 2’s , and increases with 1. The configuration 
which corresponds to (6) will persist and produce a division under 
proper conditions. At the early stages of the whole process the terms 
which correspond to damped periodical oscillations will be of the same 
order as the term (6). Each of those terms will produce temporarily 
periodical deformations of the cell. Since the periods are in general. 
incommeasurable, the net result will be a somewhat random-like pul- 
sations of the cell. These pulsations will die out and the division be- 
gin. But this is exactly what we see on micro-motion pictures of di- 
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viding cells when the process is speeded up about one hundred times. 
(C.f. the Canti film.) 

It may happen that » will also be real and positive for a value 
of 4 which corresponds to a solution of the form 


~ Iie (ar) P; (cos 0). (20) 
Vr 


In this case a three-polar division will occur, because there will be 
three maxima of osmotic pressure. 

By investigating the range of values of Dir, Dis, Qi, Giz, Gor; 
and a2. at which this can happen it may be possible to estimate the 
probability of occurrence of triple divisions. 

Finally » may be complex with a positive real part. In that case 
the difference between the polar and equatorial pressures will oscil- 
late with increasing amplitude and the possibility of a division dur- 
ing the proper phase of one of the periods should be investigated. 

It must be remarked that even for a real positive » the concen- 
trations cannot increase indefinitely, since the smallest concentration 
must still be positive. The maximum value they can reach is deter- 
mined by the value of the time-independent solution on which the 
variable solutions are superimposed (loc. cit.). Should a division oc- 
cur at all, it must occur before these limiting values are reached. This 
is likely to introduce some additional inequalities as a condition of 
division. 

The important thing in the present case is that the division of 
a cell is caused by ordinary osmotic pressure, and not by the diffusion 
drag forces. 


LITERATURE 
Rashevsky, N. 1948a. “On Periodicities in Metabolizing Systems.” Bull. Math. 
Biophysics, 10, 159-174. 
Rashevsky, N. 1948b. Mathematical Biophysics. Second Edition. Chicago: Uni- 


versity of Chicago Press. ; 
Rashevsky, N. 1948c. “A Note on Biological Periodicities.” Bull. Math. Bio- 


physics, 10, 201-204. 
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A NOTE ON THE DIFFUSION DRAG FORCES 


N. RASHEVSKY 


COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


_ .A formal derivation of the expression for the diffusion drag force 
is given. The expression involves the coefficient of diffusion of the solute 
through the solvent and the coefficient of diffusion which the solute would 
possess if it diffused through itself in the form of a gas, in the absence 
of the solvent. 


The theory of the diffusion drag forces presents an interesting 
peculiarity. The existence of the forces follows directly from New- 
ton’s third law without any additional assumptions. Since the solvent 
offers a resistance to the flow of the solute, the latter must exert an 
equal and opposite force on the solvent. It would seem that under 
those conditions the force should be easily expressed in terms of the 
diffusion coefficient which is a measure of the force of resistance of 
the solvent. Yet apparently rather complex kinetic considerations 
have been hitherto needed (Rashevsky, 1948) in order to arrive only 
at an estimate of the order of magnitude of the forces. 

The purpose of this note is to outline a simple derivation which 
must be expected to be possible. 

Let D, be the self-diffusing coefficient of the solute, considered 
as an ideal gas, in the absence of the solvent. Then, in the absence 
of the solvent, a gradient of concentration c would produce a flow 


per cm? 


Q=—D, gradc. (1) 
The mass velocity, v) , corresponding to this flow, is given by 
Q= VC, (2) 
or 
i ne (3) 
¢C 


Equations (1) and (3) give 
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D 
Vo = — — grad c. (4) 
c 


The gradient of concentrations in a gas produces an aerodynamical 
force F per unit volume equal to — grad p, where p is the pressure. 


Therefore 
F ie d (5) 
= — —— orad ¢. 
M of 


In the absence of inertial forces and in a steady state, the force 
F is equal to a frictional force produced by the velocity v, and equal 
to av), where a is a coefficient. Hence 


RT 
Woe eB, (b (6) 
or 
es ae C. (7) 
aM 
Comparison of (4) and (7) gives 
Dek 
Treeay (8) 
or 
Pelviese 
aati aik (9) 


In the presence of the solvent, the latter exerts an additional 
force f, of resistance upon the diffusion flow and makes the latter 
smaller than given by (1). The flow now is 


—Degradc, 
where D << D,. ae 


The corresponding mass velocity is now equal to 
ee 
Ue shad C3. ela (11) 


We may regard the effect of the solvent as equivalent to a force f,, 
directed opposite to F, in other words, a force of resistance added to 
av. Then for a steady state we have 


F—f,=av. (12) 
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Introducing (5), (9), and (11) into (12) we find 


kT RT ¢ D 
ee ea (13) 
or 
fe Ay a) 
= 57 (p,— 1) arate ie 


f= (1-—) grade, (15) 


Se Lae 
The quantity D. is about 10-* — 10°. Hence almost all of the osmotic 
: 


RT 
force — va grad c is sustained by the solvent. Inside of the cell we 


have D;; outside of the cell, D. > D;. If D; = D., then the potential 


of the force f , 
ali ( 1 a ) C (16) 
Ya Dea ae 


would be continuous throughout the cell and the whole space if the 
permeability h is infinite and c is continuous at the surface. Applica- 
tion of Betti’s equation to the inside and the outside of the cell would 
give zero deformation for this case. If, however, D; $ D., then 


RT ( D; ) 
- 1 —— Cc (17) 
¢ Mu D, 
and 
Rr ( DD: ) 
{= 1— C, 18) 
$ M D; \ 
and at the surface there is a discontinuity 
RTD D; 
0 Se ae C, 19) 
OD Gi beGs Te D, ( 


which permits a resultant elongation. 

Compared to previously derived expressions (Rashevsky, 1948) 
the value of the effective potential, and therefore of the effective force, 
is reduced by the factor (D. — D;) /D. ~ 10+. H. D. Landahl’s (1942) 
work indicates, however, that in Arbacia eggs the forces needed to 
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produce division are about 10* — 10-° times smaller than those com- 
puted from plausible values of the effective metabolism and from the 
assumptions that the force per unit volume is of the order of 
— (RT/M) gradc. 

Inside of the cell the potential of the diffusion force is 


gah D; 
M Ds 
outside of the cell it is equal to 
RT D-. 
b= == (1— Je. (21) 
M D, 


Let h = wo, so that c is continuous at the surface. Does the discon- 
tinuity 

D.—D; 
See SY 


D, (22) 


Ps Paes 
of the potential introduce any additional surface forces? 

To answer this question let us consider the realistic case in which 
the boundary between two phases is not a mathematical surface but 
a transitory layer of the thickness 6 in which all properties vary con- 
tinuously from one phase to the other. 

The potential ¢ will also vary continuously from ¢; to ¢-. The 
force per unit area is equal to 


ue f [1-— |Ga= 


De dx 
(23) 
6 fe) 
Ue seat Dinca 
- — daz + — a) — d 
: Ee D, } des x, 


if « is chosen in the direction of the normal. 

The first term on the right-hand side is equal to — [c(6) — ¢c(0)] 
= 0, because of the continuity of c. In the second term D(de/dx) = 
const. = Q. Hence the second term equals Q5/D, and tends to zero 
with 6. A discontinuity of D thus does not introduce any additional 
surface forces. 

Tt should be noted that expression (15) represents the force on 
the unit volume of the solvent as a whole. Any suspended particles 
are to be considered parts of the solvent. The addition of such par- 
ticles affects D , and thus affects the force. Herein is an essential dif- 
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ference between the derivation suggested here and previous ones in 
which the solvent had to be broken somewhat artificially into solvent 
proper and suspended particles (Rashevsky, 1948). 


LITERATURE 
Landahl, H. D. 1942. “A Mathematical Analysis of Elongation and Constriction 
in Cell Division.” Bull. Math. Biophysics, 4, 45-62. 
Rashevsky, N. 1948. Mathematical Biophysics. Second Edition. Chicago: Uni- 
versity of Chicago Press. 
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NOTE ON A CASE OF NONLINEAR DIFFUSION 


N. RASHEVSKY 


COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


When the function Q in the equation V2e + Q(c) = 0 is positive 
and is of a specified kind, the equation admits of a centrally spherical 
solution such that c is positive everywhere, tending to zero at infinity 
and dc/dr = 0 atr = 0. Physically this corresponds to a local concen- 
tration of the solute in an infinite medium without any membranes pres- 
ent. This result would indicate the possibility of the formation of spon- 
taneous concentrations and non-uniformities in non-linear diffusion fields. 
Possible biological implications are mentioned. 


Consider the general diffusion equation 
DV?e+Q, (1) 


where D is the diffusion coefficient, ¢ the concentration, and Q a func- 
tion of c with the following properties: 


1) Q(c) is continuous and positive for c > 0. 

2) Q(0) =0. 

a) CO Co) — Q, (| const). 

4) In the interval 0, c,, where ec, is a constant, Q(c) varies 
as ac”, where n > 8, and a is a coefficient. 

5) Inasmall interval c, — c., where c, is another constant such 
that ec. — c, << ¢,, Q(c) raises very suddenly, then again 
begins to vary slowly, and for c > ¢, tends asymptotically 
to Q, so that Q, —c << Q, fore > @. 

6) The coefficient a is sufficiently small so that ac,” << Q). 


Putting D = 1 in equation (1) without loss of generality and 
writing it in polar coordinates for the case of spherical symmetry, 
we obtain 


dc 2 de 
a gue, TO MeY. (2) 

dr? r dr 
Let us construct a solution of (2) such that when r = 0, then 
C=C, > c,anddc/dr=0. (3) 


. 15 
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Since for c > ¢,, Q(c) is very nearly constant and equal to Qo, 
the condition (3) shows that at r= 0, d?c/dr? < 0. If Q(c) were 
exactly constant and equal to Q, , then we would have c = ¢ — Qor?/65 
de/dr = — Qor/3; @c/dr? = — Q./3. As Q(c) actually slightly devi- 
ates from a constant, the actual values of de/dr and d2c/dr? will be 
slightly different, but d?c/dr* will remain negative for r = 0. 

But if at r = 0, d2c/dr? < 0, then, as r increases, both d?¢/dr? 
and de/dr remain negative, and therefore ¢ decreases. As long as ¢ 
remains above c,, Q(c) remains appreciably constant, while the abso- 
lute value of (1/7) (de/dr) increases. When however ¢ drops below 
c., then any further small decrement of c causes Q(c) to drop sharp- 
ly according to 5) and 6). Now the term Q(c) in (2) becomes neg- 
ligible as compared with dc/dr, and c varies according to the equa- 
tion 


Pe 2 de 0 (4) 
dr? “+ dr 
which, integrated, gives 
1 
C2 Corts (5) 
7 


Now we have (1/r) (dc/dr) « 1/7?, and therefore because of 
4), Q(c) will remain negligible as compared to (1/r) (de/dr), if it 
has become negligible at any one point. 

Our solution represents a more or less steep concentration of a 
solute in an indefinitely extended medium. The function Q(c) cor- 
responds to a production of the solute. Since equation (2) is an 
equation for a stationary state, we see that for a function Q(c) as 
specified above, such concentrations of the solute may exist in a sta- 
tionary state. 

We shall not investigate here whether the stationary state de- 
scribed in this paper is stable. We shall, however, point out that with 
linear diffusion we cannot obtain such solutions, stable or not. For 
Q = q (= const) the only solution which is finite at r = 0 is 


qr? 
a One eat (6) 


This does not vanish at infinity, becomes negative for r > \/6c./q 
and thus can have no physical meaning in the whole space. For 
Q(c) = k*c we have 
sin kr 
c=A ’ (7) 


ts 
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This vanishes at infinity but in an infinite number of regions is nega- 
tive and therefore again physically meaningless. Only within limited 
space, bounded by proper membranes, can solutions (6) and (7) ac- 
quire physical meaning. The solution of the nonlinear equation (2) 
is physically meaningful in the whole space and represents a local- 
ized accumulation of the dissolved substance, which does not need 
any membranes for its maintenance. 

For the case of spherical symmetry, when the concentration de- 
pends only on one coordinate, 7, we can always find a form of Q(c) 
in (2) such that the solution of (2) will be a prescribed function of 
r. For if c is given, so that 


cei (Ty, (8) 
then V2c is also given as a function of 7, 
Ae = IEG a (9) 


Eliminating rv from (8) and (9) we find V?c = v(c), and we put 
Q(c) = v(c). 

Should the possibility of such accumulations occur in other non- 
linear cases, and should those accumulations be stable, this would of- 
fer interesting possibilities of spontaneous local accumulations of dis- 
solved substances in a cell. Those may result in interesting perma- 
nent structures. 
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Using an electrical model to represent certain features of a nerve 
fiber together with a one-factor theory of excitation, an expression is 
obtained for the velocity of propagation of a nerve impulse along a nerve 
fiber exhibiting saltatory transmission. The velocity shows a maximum 
with respect to internodal length. A critical internodal length, a critical 
radius, and a critical value for the resistance of the external medium 
are required in order to make transmission of the impulse possible. 


Expressions for the velocity of conduction of an impulse in nerve 
fibers have been obtained using various models (Rashevsky, 1933, 
1938; Rushton, 1937; Offner, Weinberg, and Young, 1940). N. Ra- 
shevsky (1938, 1948) considered the case of fibers with nodes of 
Ranvier. However the explicit value of the currents was not given 
in terms of the physical constants. It is the purpose of the present 
paper to introduce a specific network and thus find the current dis- 
tribution in terms of the physical constants of the nerve; then to dis- 
cuss the dependence of the conducticn velocity on the physical con- 
stants. 

Let a myelinated nerve fiber, whose radius, exclusive of the 
sheath, is r cm, have nodes of effective electrical width 6 spaced a 
distance a cm apart. Let y be the ratio of the internal resistance per 
unit length to the resistance per unit length of the external medium. 
Let p be the specific resistance of the internal medium and let Q be 
the transverse resistance in ohms cm? of one square centimeter of 
membrane plus any’ sheaths at the node. We shall assume that the 
myelin sheath of the internode is a perfect insulator. Let the normal 
resting potential of the nerve fiber be EF’ volts. 

In order to solve our problem it is necessary to ascertain the 
current distribution along the nerve. To accoraplish this most easily 
we represent the principal properties of the nerve by the electrical 


network shown in Figure 1. 
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Whenever the node is short compared to the internode, the cur- 
rents will probably be uniformly distributed across the cross section 
of the nerve fiber and transversely along the node. In this case, if we 
ignore the effects of capacity or inductance, the geometry of the nerve 
leads us to the network of Figure 1, except for the element repre- 
senting the zero node which will be discussed below. 


DIRECTION OF PROPAGATION : 


- pa 
Eada) 


eh RS GED o 8 ee  ton ee ae Leee -- OUTSIDE 
‘ “7 oe E z 


ie Uy.) a Ly ln+) 


-- INSIDE 
Ly Inet In nat 
FIGURE 1 

In treating the problem of conduction we shall follow N. Ra- 
shevsky (1938, chap. xx; 1948, chap. xxviii) in which the law of ex- 
citation is given by a one-factor theory (Blair, 1932; Rashevsky, 
1938, chap. xvi; 1948, chap. xxiii). According to Blair’s theory, ex- 
citation in the nerve is governed by a change from the normal in the 
amount of a factor in the neighborhood of a nerve. The rate of accu- 
mulation of this factor increases in proportion to the current applied 
to the nerve and decreases by an amount equal to a constant & times 
the change in concentration of the exciting factor. The constant k 
tells one how rapidly the factor dissipates when the current is re- 
moved. It is postulated that excitation occurs when the amount of 
the factor reaches a threshold. The smallest current which can pro- 
duce excitation is the rheobase R. We shall consider the transverse 
current through a node of Ranvier to be the current causing the ac- 
cumulation of the factor determining excitation, so that this is the 
current to be used in the equations of excitation. The current den- 
sity rheobase in amperes per square centimeter will be denoted by 
Rk = f/2ar6. We shall assume that the transverse resistance w drops 
to a small enough value at the point of excitation so that the current 
flow can be calculated by assuming its value to be zero. The poten- 
- tial at the point of excitation is taken to be zero. This accounts for 
the element representing the zero node in Figure 1. 

The assumption that the transverse resistance becomes zero at 
the point of excitation is introduced for simplicity. The properties 
we shall discuss in the present paper do not depend significantly upon 
the value to which the transverse resistance at the excited node is 
reduced. A direct consequence of this simplification in the model is 
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that there is no non-propagated response possible. Thus for the pres- 
ent we shall not distinguish between the constant k’ for the exeiia- 
tory process at a point and the constant k associated with the initia- 
tion of a propagated impulse by means of an external stimulation. 

Let an impulse travel from left to right. The node at which ex- 
citation has just occurred is taken to be the zero node. Since the re- 
sistance of this element is taken to be zero, the current distribution 
along the rerve to the right of this node is independent of the cur- 
rents on the left. Due to the asymmetry of the network, there will 
be a current distribution along the nerve to the right. Of the cur- 
rents through the nodes, that in the first is the greatest. Thus the 
excitatory factor in that node will be the first to reach the threshold 
required for excitation. In the network, excitation is represented by 
E and @ going to zero. Then the current distribution shifts one node 
to the right. The rate of this shift is the conduction velocity. 

We shall now compute the current distribution to the right along 
the nerve fiber. The net flow of current down the nerve must be zero; 
thus the current inside of an internode is equal and opposite to that 
outside (Rashevsky, 1938, chap. xix; 1948, chap. xxvii) or, from Fig- 
ure 1, i, =i,’. From Kirchoft’s laws we get the equations from which 
to obtain a single equation determining the current distribution (c.f. 
Rashevsky, 1938; 1948). The junction equations give 


Lys = Ass ws ae a Dans oe. tot ’ (1) 
and the equation for the potentials may be written 
0=t(yS + S) + theta. (2) 


Eliminating the transverse currents we obtain the difference equa- 
tion 
(y +1)S + 2 


@ 


tn + Ina = 0, (3) 


Unt 


for which the solution is 
i, =A(1+ p— V2y + y*)", (4) 
where A is a constant determined by the boundary conditions and 
where | 
(1 + y) apd = (1+ y)S 
eeey rQy ST oe 
In order to evaluate the constant A we consider the currents at 
the excited region. We find, corresponding to equations (1) and (2), 
the equations 


(5) 
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L=U+h; (6) 
B=ySt, + Sia ho. (7) 
Eliminating 7,, 7,, and i from equation (7) by equations (4) 

and (6), we obtain an equation which, solved for A , gives 
A=E/of1lt+typ—V2ypt yily + Vay + y?] . (8) 
From equation (1), 4 = 7% — i.. Introducing (8) into (4) we 

can obtain 7, and 7. In this way we find 

1,=E(V2y + y?—y)/o(V2y ty +). (9) 


Using essentially the same argument as was used by N. Rashevsky 
(1938, chap. xx; 1948, chap. xxviii) we obtain an expression for the 
time required for the excitation to move from one node to the next: 


At == toe. ( ues 2ar6 (1 +papi—Fy/ Zinide mF) (10) 
i — 2nrdR 


Hence the velocity is a/At. 
We define 2 as 


E 
re (11) 


QR 


The quantity 4 is customarily referred to as the “safety factor” and 
is a measure of the ratio of the current that flows through the nodes 
of the functioning nerve to the minimum current that can cause exci- 
tation. Upon elimination of 2, from (10), using (9), and making 
appropriate re-arrangements we obtain the final result: 


2 
AVI + 27rQy/apd C1 + vy) —A—2 


v=ak/log. (: + ) eAGh 

Discussion. 

We now consider the effects of the various parameters on the 
velocity of conduction. The velocity is directly proportional to k. The 
distance between nodes, a, enters in two places, one tending to increase 
the velocity and the other tending to decrease the velocity. An in- 
crease in A causes an increase in the velocity. The parameters FE , 2, 
and F determine 1. The other factors enter as a combined quantity 
which, however, also contains 2. 

The accompanying figures show the dependence of v on several 
of the physicochemical parameters of the nerve fiber and its environ- 
ment. Unless otherwise indicated, the parameters used in making the 
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plots are as follows: k = 2000 sec, a= 0.1 em, r= 5 X 10-- cm, 
2 = 10* ohm cm?, y= 1, p= 10? ohm cm, 6 = 2 X 10-3 em. From these 
values we may calculate w and S from # = Q/2ar6 and S = p&/yar?. 
If we know the chronaxie time 7; for a given nerve fiber, the constant 
k is obtained from the relation kr = log.? (Rashevsky, 1938, chap. 
xvi). The value chosen is representative of those determined by vari- 
ous authors. To find a value for 2 we assumed that the myelin sheath 
is superimposed over a membrane whose resistance for a square cen- 
timeter is of the same order as that of the crab nerve (Hodgkin, 1947). 
The impedance of the sheath is taken to be essentially infinite (Tasaki, 
1939). Using this value of 2 and setting # — 0.05 volts, then for 
A=2, R = 1.7 X 10° amperes/cm? or 1 X 10 amperes per node; 
using only the electrical parameters of the nerve, we find, for com- 
parison, 2% — 2.3 X 107° and 7, = 3 X 10-* amperes. 

If we consider the limiting value of v when a > 6 and 6 > 0 
we find an expression for. the velocity of the non-medulated nerve 
fiber which we may write in the form 


TQ 
3 ieee (13) 
2p(1 + y) 
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This expression would hold for a continuously medulated nerve 
fiber if @ were taken to be the transverse impedance of one square 
centimeter of the membrane together with the myelin sheath. Equa- 
tion (13) is the same expression, except for differences in notation, 
as that obtained by N. Rashevsky for a nerve with a continuous sheath 
(Rashevsky, 1938, chap. xix; 1948, chap. xxvii). The same relation 
was also obtained for a continuously medulated fiber considered as a 
special case of a fiber with an interrupted myelin sheath (Rashevsky, 
1988, chap. xx; 1948, chap. xxviii). The above expression is also for- 
mally the same as that given by W. Rushton (1937, cf. Offner, Wein- 
berg, and Young, 1940). 

If we obtain v for a > 6, the value is somewhat smaller (v = 4.4 
m/sec instead of 4.5 m/sec) than the value for 6 = 0. The differ- 
ence arises because of the implied interference with the current flow 
of a zero length internode. 

In Figure 2 we consider the variation of v with the internodal 
distance, a, for three values of the safety factor 4. It is seen that at 
a—6,v has a value determined by J; it increases to a maximum as @ 
increases, and finally falls to 0 for a critical value of a. It it were 
possible to effectively insulate a node, such behavior could be investi- 
gated experimentally. Since the critical value of a depends upon jd, 
this technique might be used to evaluate 4 for various nerve fibers. of 
the myelinated type, or at least impose a restriction on its value. 


VELOCITY IN M/SEC. 
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In Figure 3 we plot the value of v as a function of the nerve 
radius, The radius considered should be the radius of the cylinder 
occupied by the axoplasm. In general the velocity increases with r. 
This fact has been reported by many investigators (Erlanger and 
Gasser, 1937). However it is seen that a critical radius exists, i.e., 
for a certain set of physical parameters, the fiber must be greater 
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than a determined radius for propagation to occur. For q = 0.1 cm 
this critical radius is seen to be 0.05u for 4 = 2 and 0.16u for 41= 1 
For a= 1 cm, 4 = 1 the critical radius is 1.5u. 


VELOCITY IN M/SEC. 
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In Figure 4 we consider the effect on the conduction velocity of 
the electrical impedance of the medium in which the nerve is im- 
mersed. This is a factor which may be varied by placing the nerve 
in an oil medium, by adding a non-electrolyte such as sucrose to the 
external medium, or by imposing physical constraint on the flow of 
the external current. It is seen that as y decreases, that is, as the 
external resistance is increased, we have a smaller conduction veloc- 
ity. The interesting feature of this curve is that for the non-myeli- 
nated fibers there is no value of y for which v = 0. However, if the 
nerve is myelinated the theory predicts that there is a value of y for 
which nerve conduction is not possible. For example, if 4 = 1, the 


. 1 
nerve should cease conducting if placed in a medium such that y > ae 


However, if the safety factor is greater, we see that the action of the 
myelinated nerve approaches that of the non-myelinated fiber, and 
the different behaviors may be experimentally indistinguishable. 

The form of the curves suggests the following experimental tech- 
nique for the evaluation of 4 and one of the factors 2, p, or 6. If 


A> 4/(V1 + 2/y —1) (a condition which is satisfied for the typical 
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nerve in a normal environment), we may simplify the expression for 
v by expanding the logarithm. It is readily seen that the following 
relation is obtained: 


an) (ER) 
2\—+1])] =A ae 
ak apo i Ds nie 


Therefore if we know the value of k from consideration of the chro- 
naxie time, and if we measure the value of a by direct observation, 
we can find 4 and \/2rQ/apé by plotting the quantity 2(v/ak + 1) 
against V/y/(y + 1). It must be remembered that the only values of 
y allowable are those which satisfy the inequality above. The data 
should lie on a straight line whose slope is 4\/2rQ2/apé and whose in- 
tercept on the Vy/(y + 1) axis is /272/apd. The intercept, then, al- 
lows one to find the quantity Q/pé using the measured values of r and 
a (provided one is able to make a satisfactory extrapolation to the 
axis of the experimental points). If p and 6 are estimated indepen- 
dently, then 2 can be calculated. Since 2 = E/RQ — 1, we can now 
find the value of R from a measurement of the potential EF’ of the 
nerve fiber, which in this model is the resting potential. 

In all of the above graphs the various parameters were chosen 
as representative values. If in a particular case they are indepen- 
dently measured then the velocity can be computed and compared with 
the observed value. If now, for example, y is varied, but any changes 
in the other parameters are measured, then the velocity can again be 
computed. At the value at which conduction is just blocked, for ex- 


ample, substitution of the new values of the parameters must satisfy 
the relation 


2096 (1'+ y) + A) =27rQy. (15) 


Actually the left-hand and right-hand members need only be very 
roughly equal. 


Summary: 


1. A model for a conducting nerve is introduced. From this 
model the velocity of a nerve impulse may be deduced as a function 
of the physical constants of the nerve fiber and its environment. This 
expression is derived. 

2. If we consider the internodal distance of a myelinated nerve 
as the only parameter of the system, it is found that there is a maxi- 
mum velocity and that there is a critical distance between nodes. For 
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internodal distances greater than this critical distance, there can be 
no propagation. 

3. Similarly, there is a critical radius for a myelinated nerve 
fiber. Considering only the variation of the radius, it is found that 
there exists a radius which must be exceeded by a fiber found in na- 
ture if it is to be functional. 

4, If we consider a single myelinated nerve fiber placed in a 
medium whose resistance can be varied in such a way that the prop- 
erties of the nerve fiber are not changed, there will be an impedance 
for the external medium which can prevent the nerve from conduct- 
ing impulses. 
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The steady state kinetics of some typical catalytic systems of biologi- 
cal importance have been formulated. The conditions for the existence of 
a maximum or limiting velocity are examined and discussed. In particu- 
lar it is shown that the limiting velocity for a given component is simply 
the rate expression for a given number of steps of the overall process; 
from the general condition for a limiting velocity these steps may be 
specified. The stringency of the conditions which must be imposed upon 
the steady state solution in order that it may be assumed that one or 
more steps are essentially at equilibrium is pointed out. The application 
of the general method to coupled or branched systems and to cyclic sys- 
tems is briefly discussed. 


Many biological processes consist of consecutive reactions. Wheth- 
er the process in question is a single enzyme system which catalyzes 
the oxidation of a substrate through a so-called “hydrogen transport 
chain,” or the removal of the products of reaction from a complex 
system by another system or subsequent sets of systems, this feature 
of consecutive rate processes makes its appearance. Most frequently 
such reaction nets will be found in the steady state because various 
intermediates reach a stationary or time independent level and do not 
accumulate indefinitely. In general this time independent state will 
be different, and often far removed, from the equilibrium state. In 
- extreme instances the conditions which specify the steady state may 
approximate the equilibrium condition for some of the steps in the 
overall process. We exclude here the case in which all steps are in 
equilibrium which would require that no net change occurs in the sys- 
tem. The prominence of consecutive processes in biological systems, 
along with the fact that the solution of the generalized set of consecu- 
tive reactions is not available unless the conditions of time indepen- 
dence may be assumed, makes the analysis of steady state kinetics a 
coordinate feature of many biological problems. Furthermore, in view 
of the fact that in general it is not valid to assume a priori that any 
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of the steps in the process are near equilibrium, it is of importance to 
initially examine true steady state solutions. 

The classical Michaelis-Menton formulation (1913) of enzyme 
catalysis rests upon an assumption which, as will be pointed out, may 
seldom be valid. When this assumption is removed, the formulation 
applies at best to a simple two-step process although it has been 
widely utilized analogically for the description of multi-step systems. 
The ‘“Michaelis-Menton constant,” K,, a perfectly definite concept 
in the simple theory, has often been taken to have the same physical 
significance in multi-step systems. The kinetic description of complex 
systems is frequently formally identical to that of the simple systems 
because there is a limiting velocity and an “apparent K,,”’ which char- 
acterizes a hyperbolic relation between rate and substrate concentra- 
ticn. It is clearly desirable to examine the make-up of the constant, 
which makes its appearance as the counterpart of K, in multi-step 
systems, and to inquire into the conditions for a limiting velocity. 
Where it is feasible, the stringency of the conditions which must be 
imposed upon the steady state solutions in order that the assump- 
tion of equilibrium in some of the steps may be justified, will be ex- 
amined. 


The General Formulation. 

The formulation which at once recommends itself is that given 
by J. H. Christiansen (1935) and generalized by L. P. Hammett 
(1937; see also references in Hammett). This formulation has been 
used by M. F. Morales (1947) in discussing limiting reactions and 
has been applied to the blood coagulation system (Hearon, 1948). It 
is still remarkable that this method of choice was not applied earlier 
and more extensively in biology. For definiteness and as a basis for 
reference in what follows and in, forthcoming papers, some details 
of the method will be given. 

For the general set of reactions 
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Ga Gs —X,+B,, 
C.+X, =X.+B,, 


G,+Xi4 =X + B;, 


Ca Ae as =—B,, ’ 


denote by k; and k_; the forward and reverse rate constants for the 


ki 
ith step and let K; = —. If by definition 
w,=k, [C,] [C.], 
w=k, [Ci]; +471, 


(1) 
wi=—k, (Bil; t1# 18, 
W_n — iD ss > 
eee 0 
then the steady state conditions, ——————, lead to 
PS Ww ee} 
1 = Wo =e W_2Ke > 
(2) 


1 = Widia — W-i0; , 


i Wyn — W_nOo » 


=-,o= EX] , and » is the overall rate of conversion of 


where o = 
Vv 
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C, and C, to B,. It must here be assumed that the concentrations of 
the reactants C; and products B; are maintained constant. This con- 
dition is usually met by having the concentrations of these constitu- 
ents sufficiently large so that their change during the time of obser- 
vation is inappreciable, or in other ways as has been discussed by A. 
C. Burton (1939). There is still another way particularly pertinent 
to biological problems in which this condition of constancy may be 
satisfied: If a reactant in one reaction is regenerated in a subsequent 
reaction, the substance is conserved to the system and its concentra- 
tion maintained stationary. If the last step is irreversible, this re- 
striction of constancy may obviously be removed from B,. The solu- 
tion of equation (2) is 


n n 
IT Ua II w-.; 


v 


(3) 


WoW, +++ Wy + WaW3 +++ Wy + .506 + W_4W_2 ++++ W_(n-1) 


From equation (3) it is seen that if any w_; = 0 the second term in 
the numerator of v vanishes. In what follows, we shall often assume 
that w_, = k., = 0. The extension of the solution (3) to biological 
problems involves a feature of difficulty not always encountered in 
chemical kinetics; viz., the application of the material balance to 
(catalytic) components which are conserved in the system. If C; is 
such a component, then 


[Ci], =([Ci] +> LO], (4) 


where [Ci], is the total concentration of C; , and X;°: is any interme- 
diate which contains C; in bound form. Condition (4) may be im- 
posed on the problem by including this equation in the set (2). In par- 
ticular cases this is the simpler procedure. Alternatively, equation 
(2) may be solved for [X;°:]; this expression, which in general con- 
tains v, substituted in equation (4) gives [C;] as a function of the 
ws and [C;i],. The solution of equation (2) for any [X;], first given 
in equivalent form by M. F. Morales (1947), may be written 


+ w 1 t, 
[XJ=0—*f1—9 | —+ WwW 
kK W_k W, WyWs 
5 
ie eee | . 
WW +++ W; 


If the ith step is separated from the 1st step by a number of steps, 
equation (5) and its use in (4) are not simple; if the 7th step is sepa- 


JOHN Z. HEARON 33 


rated from the mth step by a smaller number of steps, a shorter form 
is available. This expression, which will not be derived here, is 


[xg=m {1+ 0 : a) gaeecases 
=n Wk Wy, W-(n-1)W-n 
W (i+2) oH Wy | (6) 
+S 
Wea ab. | 


Equation (6) is the result obtained from equation (2) by determinants, 
or by solving the nth member of (2) for an , substituting this into the 
(n—1) st member, etc. It should be noted that the quantities in brackets 
of equations (5) and (6) are such that they may be written at once 
from comparison to (3). In equation (5) the expression in brackets is 


1 

simply V0” where V,,; is the rate of steps 1 through 7 terminating 
1,7 

irreversibly in i. Similarly, the expression in brackets of equation 


1 
(6) is ———, where V. -n (41) is the reverse rate of steps n through 
—n , (141) 
4 + 1 terminating irreversibly in — (7 + 1). On this basis, Vi; and 
Vin Gay may be readily constructed, by using equation (3) , in more 
convenient form. In the event that k_, = 0, equation (6) becomes 


Wize 20° Wn i W-(i41) Wiss 2° Wn TF isees W_(i+2) °°°* W-n(n-1) 
[Xi] =o (7) 
Wii tts Wn 
= ON A Fre ’ 


where Gi: is the function such that [Xi] Gian = Vian and is ob- 
tainable from equation (3) by setting up Vian and deleting [Xi]. 
These properties will be convenient in discussing the change in [X;] 
introduced by changing parameters of the system. In general equa- 
tion (5) will give i terms in the bracket; equation (6) will give n—1 
terms. The choice is obvious in a specific case. 


The Simple Enzyme System. 
It will be convenient as a basis of comparison to solve the simple 


_system . 
E+S=2E-S; 

Bes 2b. 4 ps 

where E, S, and p denote enzyme, substrate, and product respective- 
ly. The solution is 
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WW. _ hikes [EF] [S] (8) 

W4 + We Ky + hea 

using equations (7) and (4) [Eo] = [EZ] + [E-S] 


v 
[i -S)=—; (9) 
We 
Ko ae Dg i [es 
y= He Cod LS) _ fe Wa) 181, on 
is = ie i ae LS] 
oes ee 
in he SPF : : 
Here corresponds to K,, and it is readily seen that if 
Keg : . 
2 Sey A ea rag = ra which corresponds to the assumption that 


the first step is in equilibrium. When the condition k.. >> kz is satis- 
fied, then the determination of K,, from rate curves gives the recipro- 
cal of the equilibrium constant of the first step. Under these condi- 
tions 1/K,, is a measure of the affinity of the enzyme for the substrate. 
It is noted that the above is the exact equivalent of the Briggs-Haldane 
formulation of the problem (Haldane, 1930). This case is given here 
for its illustrative value and for comparison to complex systems. 
While it is not unreasonable that there should be cases in which 
k, << k., it should be noted that in the few instances where indi- 
vidual determinations of k,, kz and k., have been made, the results 
have been quite the converse; i.e., kK. << k, (Chance, 1943; Gold- 
stein, 1944). In such cases the constant, which would be observed as 
Ky , would not be an equilibrium constant. In the event that k, >> ky 
and K,, ~ k./k, , the constant K,, would exhibit a temperature depen- 
dency typical of an equilibrium constant; or, if AF.* ~ AF,t, would 
exhibit practically no temperature dependency. If [S] > ow, from 
equation (1); v > k,[E,]. Under these conditions [EZ] — 0 and 
[E'-S] > [E,]. Thus the limiting velocity is the rate of breakdown 
of [E'- S] via step two; this is said to be the limiting step and [E] is 
said to be the limiting component. In this case, the hyperbolic relation 
(10) between [S] and v is clearly due to the saturation of the enzyme, 
and is a direct consequence, analytically, of imposing the material bal- 
ance (9), because when [E,] >> [E-S], v is given by equation (8) 
with [E'] = [E,]. 

That the existence of a limiting velocity is not necessarily due 
to the saturation or binding of a given component is seen as follows: 

For a system of n= 3, equation (3) becomes, if w.,=0, 
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WiW2W2 
er ; (41) 
W2Ws = W34W3 =F W_4W_» 


If no material balance is necessary, equation (11) gives the correct 
steady state velocity, and we are assuming that no C; is in the bound 
form analogous to E-S. From (11) it is seen, recalling the defini- 
tion of w;, that v is hyperbolic in [C.] and [C;]. The limiting ve- 
locity is, when [C,;] > o, 
W1We 
Cs ey, (12) 
We + Ws 

which is recognized as the velocity of the first two steps with w.—=0. 
The interpretation is clear: The value of [C,] can effect the overall 
rate of the system only by causing the consumption of X. and thus 
repressing the reversal of step two. From equation (7), when 
[C3] > 0, [X.] > 0. The limiting velocity is the rate at which 
steps one and two can supply X, to step three, w_. being zero. This 
interpretation is general for such a system and is intimately related 
to the fact (also pointed out by Burton, 1939) that if any step is 
made irreversible, the rate becomes independent of the following 
steps. Thus if the concentration of the kth reactant is increased in- 
definitely, the rate becomes that of the first k — 1 steps. In the case 
being discussed, no mechanism has been provided for a saturation ef- 
fect and the limiting velocity cannot be due to a limiting component. 
Furthermore, in general, no single step becomes limiting; however, 
if [C.] ~ o,v becomes the rate of the first step. It will be noted 
that v is a linear function of [C,] and no limiting velocity exists for 
C,. It is clear from equation (11) that if [C.] alone is varied, the 
quantity which would appear as K,, is 


Ws z W4W-2 . 
kp kW 
and similarly for [C,] the apparent K,, is 
W4W» 


kz (We + Wa) 


If a component in the system is conserved, the above discussion 
must be modified. In particular, if a reactant in the first step is re- 
generated in some subsequent step, its concentration will appear in 
the w for the reverse of that step. It will then be true, as is evident 
from equation (3), that a limiting velocity exists for C, and Co. Fur- 
thermore, a material balance (4) is necessary for a conserved com- 
ponent. These points will be covered by treating a concrete example. 
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A Typical Dehydrogenase System. 
Consider the system 


E+S=E£z-S 
E-St+Co@E+Co-H+p 
P+Co-H—~P-H+Co 
for which 
w,=k, [E] [8S] waka 
We he Coq W-2 =k» [E] [p] (13) 
w, =k; [P] ; W-3=0; 


where Co and P are coenzyme and flavoprotein respectively and Co- H 
and P-H reduced coenzyme and flavoprotein respectively. In general 
such a system operates by having P- H reoxidized by other “carriers” 
which are ultimately reoxidized by O.. This case has been solved, 
but those steps introduce no new features and are omitted here. The 
complex -S-Co has been omitted because its inclusion introduces 
only additional rate constants. The equations as written amount to 
the assumption that reduced coenzyme and oxidized substrate leave 
the EH — surface simultaneously. The expression for v is 


ies Iktok, [EZ] [S$] [Co] [P] 
~ Keekts [Co] [P] + kak, [P] + kak [E] [1° 


Equation (14) is the correct expression provided [E'] , the concen- 
tration of free enzyme, is inserted from 


[Bok = LE] Oo Lhe Ss] 


(14) 


‘ (15) 
[EZ -S] =n LED ES I. es 


for from equation (15), [EZ] is dependent upon all other concentra- 
tion terms except when [S] is small. If (14) is written as 


A [E] 


"~~ BoC tee ee 
where 
A=k,k.k, [S] [Co] [P], 
B=k; [P] (k, [Co] + k.), ety) 
C=k.k. [pl], 
then 


+Strictly speaking, it would be necessary to introduce the possibility of the 
formation of E:Co , E-Co-H and E-p-Co-H with Co, Co-H and S Ae p Eharanters 
ized by different affinities for the H — surface. If S and Co react sufficiently fast 
on the H — surface, the concentration of E-Co may be ignored. 
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ste =1 [Bb [S] [Z] + Ch, [S] [E}?—A a A 
len B+C[E] LPO: 
and combining the two members of equation (15) gives 
A 
(=~ 2x, [SP-C [hel —B) 
[Zz] =— 
2G (1 FG-[S]) 
(19) 


A 2 
\ (x. [S] Be igen [2.1 + 4B [E,] C(1 + K, [S]) 


+ 
2C (1+ K, [S]) 

When equation (19) is used in (14) it is seen that v is non-linear in 
[S] and a limiting velocity will exist when [S] ~ o. Perhaps more 
important it should be noted that v is not directly proportional to 
[EZ] , the total enzyme concentration. The non-linearity of equation 
(14) in [EF] is clearly due to the fact that EH is regenerated in the 
second step and the non-linearity of [S] results from the imposed 
condition (15) which takes into account that EH is bound by S in the 
form E-S. The fact that v is proportional to [E,] , a property of 
the simple system expressed by equation (10), is often assumed to 
be true for complex systems. That v is not directly proportional to 
[E.] is a common finding (e.g. Hearon, 1944; Dewan and Green, 
1938). As will be evident from the system being discussed, the depend- 
ency of v upon [E>] will be conditioned by the particular steps in which 
E enters and is regenerated in the system. From equation (19), when 
[S] — o, [EZ] — 0; we will need, in order to examine the limiting_ 
velocity in this case, the expression for Lim ([E] [S]). By fac- 


[S]00 


toring 
A 
——— — Bo 
(ex, [S] i= +B-—Cl[ 1) 


from the radicand and expanding the result, the value of [E] for 
large [S] is found to be 


eB el 
a TS] {BKi— 2}, 

(20) 
en 5) 


a1 


The limit sought is thus 
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Lim ([E] [S]) = BLE.) 23 [Ee] 

ee SBR, pater 
Substituting equation (21) into (16), with [ZH] — 0, gives as the 
limiting velocity when [S] ~ 


Wee) = key [Co] [Eo] a (22) 


Equation (22) gives v when [E'-S] =H, and w.,=—0. This is readily 
perceived physically since when [S] ~ o, [E'-S] > [EZ] requiring 
[EZ] — 0, the second step is rendered irreversible. In an entirely 
similar manner it may be shown that when [Co] > o, [EH] ~ [Eo] 
and the limiting velocity is 


Vn =k, [S] [Fol , (23) 


which is the forward rate of the first step when [ZH] = [EF]. This 
is seen to result physically from the fact that when [Co] 7 o, 
[E-S] — 0 rendering the first step irreversible; the overall rate of 
v is then independent of the subsequent steps. When [P] ~ o it can 
be shown that [H’] approaches 

[E] = (k, [Co] +k.) [Eo] 

: k, [Co] +k, [S] +k. 
Substitution of equation (24) into (14), which for large [P] has the 
form (25) ; 


(21) 


(24) 


_ Take [E] [S] [Co] 


| COT Sea) ea ae ) (25) 
gives 
aioe 
v.00) — IE IS Teel, (26) 
re + [S] oe aa [Co] 


which can be shown to be the rate of the first two steps with w..=0. 

The results (22), (23), and (26) show the extent to which the 
principles discussed with respect to equation (11) apply to a system 
in which one or more components are conserved and the modifications 
which are necessary. The results may be generalized as follows: There 
will be a V,, for a component occurring in the first step provided that 
it (1) is conserved in the system and regenerated in a subsequent 
reversible step or (2) is capable of binding the free form of a cataly- 
tic component. For any C;,i # 1,710, which does not bind a cata- 
lytic component the corresponding V,,‘°) is the rate of the preceding 
1 — 1 steps as already noted. If any C; does bind a catalytic compo- 
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nent, the corresponding V,,‘° is the rate of the subsequent steps 
a+ 1 through i + &, terminating irreversibly in the (4 + &)th step, 
where the free catalyst is regenerated in the (2 + &)th step. In the 
rate expression for the (¢ + 1)th through (i + €)th steps, the con- 
centration of the bound form’ of the catalyst in question must obvi- 
ously be set equal to the total concentration of the catalyst. These 
general principles, which are difficult to state verbally, may be imple- 
mented with the following extension of the case already treated: 


(3) tel oS Sh = Ss 

(b) CotE-SE-S-Co, 

(c) E-S-Co@E+p0+4+0Co-H, 
(dyes Co-e =P P.7-+ Co. 


Here FE -S-Co has been included and its concentration will be regard- 
ed as appreciable relative to [H,]. The quantity V,,“ will be the 
rate of steps (b) and (c) with [H.-S] = [E,] and w., = 0, since E 
is regenerated in step (c). The quantity V,,(°” will be the forward 
rate of step (c) with [E-S-Co] = [E,]. If additional steps were 
interposed between (c) and (d), these steps would be included in 
Vn” down through that step regenerating #. The quantity V,,” 
is the rate of steps (a), (b), and (c) with w., = 0. This example 
includes all of the cases discussed generally above. 

These deductions are of importance for the experimental analysis 
of such systems because they state that the various limiting velocities 
are given by the rate expressions for a discrete number of steps, the 
steps are specified and their dependency on the [C;] given. Although 
in the case of the dehydrogenase systems, the individual steps or small 
number of steps may be separately analyzed, the above results indi- 
cate that it is not required that such separations be experimentally 
feasible. 

To return to the general solution (14), with [#] from equation 
(19), it must be noted that v is not hyperbolic in [S] . However, for 
large [S] , when [Ff]? << [E], the solution for [Z] is 

(El = pe ie) By nore ; (27) 
; B+ kyks [P] ES C [Eo] 

and equation (27) put into (14) gives v as a rectangular hyperbola 
in [S].+ For still larger [S] , v is, of course, Vn®, and for small [S], 
[E] = [E.], and v is linear in [S] and [£.]; at intermediate values 

of [.S] the full expression must be used. 

B 

+This result would be formally identical with (10) but with “K,,” = ik [PI] 
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It is of some interest to inquire under what conditions it may be 
assumed that the first step, FE + S = E-S, of the system be in equi- 
librium. From equation (19) it is seen iat if C [E,] predominates 
the other termis 


fe Ves __ tel (28) 
1 aekAs) 
which is the equilibrium expression for [FE]. This is a direct mathe- 
matical result and the manner in which the physical requirements 
may be met is another question. Recalling the definition of A, B, and 
C in equation (17) and noting that 


2. — BK, [S] =— kyks [P] [S1, 


it is seen that equation (28) may be obtained on several bases. For 
example, if k, is large or k. and k are large, if k; is small, if [P] 
is small, etc., (28) is justified. All of the possible conditions have the 
factor in common that they amount to diminishing those quantities 
which determine the breakdown of E'-S via routes other than dis- 
sociation to H and S. Thus the assumptions which justify equation 


ul 
(28) are basically the same as those which make K,, = ae It can 
al 
also be shown that when [S] is such that [EZ] = [EZ -S]= [E,]/2, 
Vin ® 
then v # Spee: In other words, 


[E -S] v 
[Ey] vec) 
except when the values are 1 or zero, and the percentage “saturation” 


of the enzyme does not correspond to the percentage attainment of 
V,®. This is readily shown by considering that, from equation (15) 


W-4 
v=w,{1—“ w-s]| 
W, 


(29) 
el | k.  [E-S] 2 
=w,)1—- ; 
k, [SS] [E] 
and that when [E] = [E-S] = Le = [S] is found to be, from 


equation (19) 


JOHN Z. HEARON Al 


B+ = [B) 
[S] = G ; (30) 
a [EH] + kk; [P] | 


Equation (30) put into (29) gives 


> | 
(1—K,’) = (Bol + kak, (Pi + kale [P] {Co] 
exch : iy (31) 
B+ [Bu] 


The ratio of v from equation (31) to Vm from equation (22) clearly 
depends upon the values of all concentrations except [S] and every 


k;. It should be noted that when [EZ] = = [Bol and when equilibrium 


; ee 
is assumed in the simple system, the value of [S] is x From equa- 
1 


1 ; 
tion (30) it is seen that [S] = =~ only if the assumption leading to 
C >> B can be made; furthermore, from equation (29) the conditions 
E-S pe 
Les 1; Whies 
[E] 


- are incompatible if v #0. 
If in the system under discussion the “carrier” P is omitted, the 
system becomes 


Ei Sh Ss 
CotHE-S>?E+0+Co-d; 


only under the condition that during the time of observation the sys- 
tem is sufficiently far from equilibrium so that w.—=0. The solution 
is identical with V,,) already given in equation (26). This corre- 
sponds to a common experimental procedure; viz., in a system of EL, 
S, and Co, the rate of appearance of Co-H is followed optically. E. 
Negelein and E. Haas (1935) have made an analysis under these con- 
ditions of the Zwischerferment — Co;, — glucose — 6 — phosphate sys- 
tem. They state that the rate of appearance of Coy,-H follows 


v=R [Con], (32) 
and establish empirically the dependency of K upon [S] , [Cew] and 
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[E,]. Comparison of equation (32) with (26) shows that R should 
be given by Q 
pe k. [Eo] [S] (33) 


ke 
CS ae = he Hee [Corr] 


from which it follows that R is directly proportional to [Z)], a rec- 
tangular hyperbola in [S] and R~ a linear function of [Coy]. Fig- 
ures 1 and 2 show R and R~ as functions of [Coy] and F as a func- 
tion of [S] and [E,]. The correspondence found here is apparently 
satisfactory. It must be mentioned, however, that from the data it is 
by no means certain that the conditions of the steady state were fully 
met; ie., during the time of observation appreciable Coy was con- 


(s]= 4.6 
{E,j= ,0015 


[s] = 4.6 
[COn]=.05 


[Eo]=.003 
[COn]=-3 


re ee ee ee 
{Cox} OS WS ASK ZOle 25° 30 Fie Ag ier 
FIGURE 1 FIGURE 2 


verted to Co,-H. In equation (33) we are identifying [Coy] with 
the total concentration of Coy which is not strictly true in this case. 
The assumption that w.. = 0 can be fully justified from the data. It 
might also be mentioned that O. Warburg and E. Negelein (Negelein 
and Haas, 1935, p. 212) treat the data on the basis that E is in equi- 
hbrium with Coy, and Con-H and that the equilibrium constant is 
identical for the two species. The rate is then made proportional to 
[Coy] and [E,]. This analysis gives R in a form which is indepen- 
dent of [S] or, if [S] be included in the proportionality constant, is 
linear in [S] . The data in Figure 2 would seem to make this assump- 
tion inadmissable. 
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The Dilution Effect. 


If an enzyme system or tissue brei is diluted and the ratio of the 
initial and final velocities compared, rather “anamolous” results are 
obtained. Of course, in the case of tissue preparations, the elution of 
diffusible components away from components fixed in the large par- 
ticles and similar effects must be taken into account (Potter, 1941). 
Since rather unlikely explanations based on termolecular collisions 
have been given (Krebs, 1942) for the so-called dilution effect, it is 
of interest to inquire what type of dilution effect would be required 
by the rate equations being discussed. Taking equation (26) as an 
illustration, if the system is diluted to twice its original value the 
ratio of initial and final rates is given by 
i ies C bis 1 
Poe 


vy 1 


L 
5 ({S] 
v4, 


i (34) 
[Ss] + = [Co] + 1/K, 


If 


key 
if 


2 (Col <<1/Ke. -0,/0,== 8, 
ky, 
' In general, 4 < v,/v2 < 8. 

Direct and simple results cannot be obtained for the general so- 
lution (14) but it is clear that on the above basis dilution effects ap- 
proaching 32 can be obtained for a 2-fold dilution of systems involv- 
ing only four C;. The ratios obtained by Krebs using the /-amino 
acid oxidase system ranged from 2 to 19. (See for example, Green, 
1940.) The dilution effect is of further interest in more general prob- 
lems. For example, it has been suggested (Potter, 1942) that growth 
outpaces the syntheses of oxidative enzymes in tumor tissue thereby 
causing a dilution of these enzymes relative to normal tissue. Of im- 
portance also is the subject of the kinetics of the dilution effect ; i<., 
the time course of transition to a new steady state which will be re- 
served for forthcoming papers. This problem has been solved for cer- 
tain cases by A. C. Burton (1939) and discussed by him in connec- 


tion with “overshoot.” 


[S] + 
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The Turnover Number. 


In the system for which equation (14) is the solution, the mole- 
cules of Co are being “used over again.” The ratio 
v 


Coco = [Co] 


(35) 


gives the so-called turnover number. ¢,. has the dimensions t+ and 
is the number of times per unit time a molecule of Co is reduced and 
oxidized. If, for simplicity we assume that the system is operating 
at large values of [P] , then from the equation (26) 


: Ey 
fie SES EE ieee (36) 


[od celia a [Co] 


equation (36) gives fc, in its dependency on [S], [E.] and [Co].+ 
The quantity Co, exhibits a maximum value of 

ke [S] Eo 

[S] + 1/k, 


and approaches zero as [Co] — co. In general, ¢ for any catalytic 
component will depend upon the concentration of every other compo- 
nent and values of ¢ are not comparable except under strictly compar- 
able conditions (cf. Straus and Goldstein, 1943, p. 580). 


Coupled or Branched Systems. 


At+B=X, 
Ky = Xe +) eee 
Ae Een 35 0% 


For the above systems, the velocity of the »-producing system is 


kkk, [A] [B] 


v= , ot 
kok, + kik, + kk» [p] oe 
and that for the p-consuming system is 
Vi»? [p 
2= ee (38) 
K + [p] 


+In an excellent study of this kind (Corran et al, 1989), ¢ for Straub’s flavo- 
protein functioning with the lactic-dehydrogenase system has been shown to be 


a decreasing function of its own concentration as would be predicted from the 
proper modification of (14). ; 
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but the condition determining [p] is 


Veep | 
Ka [X41] — ke [X2] [p] =————, 39 
1 ELEN Neb 
from which 
ko eo Vint [p] 
Wy es U1 pe ae, y = 5" Sy 
= tae ae i (40) 


with the obvious requirement that the left-hand member of equation 
(40) does not exceed V,,“”. This requires 
kyW, 
ee ya Git (41) 
W-1 
Substituting v, from equation (37) into (40) gives for the steady 
state value of [p] 


= Vas? —4 aa; 
le scoaaes to a 
2a, 
— Ay — [A] [B] {hi kk? ae eahes a ity Fetgheny 
we WV ig = Kok, [A] [B]) 


a, = kak. (keK, [A] [B] — Vn) < 9 
dz, = K (kK, — k.?K, — kak;) [A] [B] . 


Equation (42) put into (37) or (38) gives the rate of either system 
in terms of the parameters of the other. With the proper modifica- 
tion the above system can be made to correspond to the so-called co- 
enzyme-linked systems (Green, 1940; Dewan and Green, 1937) where- 
in Co is reduced by one enzyme substrate system and reoxidized by 
another. The net result is the oxidation of the first substrate by the 
second. This corresponds to the “mutase systems” of D. E. Green et 
al (1937). One characteristic of such systems is that the rate is not 
linear in total enzyme concentration. A similar treatment can be 
made of the true mutases such as aldehyde mutase (Dixon and Lut- 
wak-Mann, 1937). In general p may enter the second system in any 
step and the solution for this system need not be the simple form 
(38), but the method of solution remains the same. 


(42) 


Systems Involving Dissociable Prosthetie Groups. 
Consider the system 
E+a2k’ 
-E*+S=2Ez-S 
foo dr eG 
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in which E is the protein moiety of the enzyme which is active only 
when associated with the group “a” to form E*. We wish to examine 
the steady state solution for this system and to inquire what condi- 
tions must be imposed in order that a study of the rate as a function 
of [a] shall yield the value of the dissociation constant of #*. Fur- 
thermore, the example will introduce a feature which obviates the 
direct use of equation (2) and its solution (38). Equation (2) will 
d (E*] 


not correctly express the condition = 0 for [E*] is reformed 


in the last step. This situation in which an X; is regenerated in a 
subsequent step has not appeared in the cases discussed previously. 
The steady state condition for [H*] is 

d [E*] 


ap UE] [a] — ka [2] — fe TB") [5] 


43 
+ Mee LE. S| he =) Oe =) 
Because of the term k; [E*-S], equation (43) is not equivalent to 
the corresponding equations which gave rise to (2). Equation (43) 
must be solved with 


LA") TS] _ (2°) bs) 


a. S|= 
[Ee S] Geet wou (44) 
Kee 
: ad[k*-S 
obtained from ae = 0, subject to the condition 
[Eo] = LE] + [2*] + LAST: (45) 
The solution for [E*-S] is 
E,| [S 
K [a] + [S] [a] + — 
Ky 
and v = k; [E*-S]. 
From (46) it is seen that if [S] << K , then v is given by 
__ ks [E.] [S] [a] 
ae eae (47) 


~ K({@] +1/K,) ’ 


and 1/K, may be determined from a study of v as a function of [a]. 
From this it appears that conditions may always be arranged such 
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that (47) is valid regardless of the inherent values of k, and k.. (cf. 
Hogness, 1942). 

If the problem is solved on a strict equilibrium basis by solving 
equations (44) and (45) with the mass action expression (48) 


[E"] =K, [E] [a], (48) 


the solution is identical with (46). This result, which means that 
the first step may be assumed to be in equilibrium, is general for this 
type of system. It is seen, in this case, to be a consequence of the 
requirement 
d [E*-S] ay; 
dt 


causing (43) to reduce to (48). If the formation of E* - S is followed 
by any number of steps, any one of which leads to the restitution of 
E* to the system, the requirement that 


d[E*-S 
Lebar 
dt 


and = 0 for any [X;] involved in these steps will reduce the 


d [Xi] 
dt 
equation corresponding to (43), to (48). For examples, the use of 
studies of v as a function of [a] to determine K, where the equiva- 
lent of (47) is assumed to hold, the reader is referred to T. R. Hog- 

ness (1942) and F. Lipmann (1939). 
Another system which must be treated in the manner here dem- 
onstrated is 
ie i, 
E,=2E.:S 
Ey: S > Eu 


where E; and E, are inactive and active forms of the enzyme. The 
solution may be obtained from equation (46) by setting [a] —1. 
Strictly, one should include the possibility that Z;-S is formed, but 
this condition is easily included in (45). More generally the above 
system would consist of a complex chain following the formation of 
E, from E; , equation (44) would then be replaced by the proper form 
of (5) or (6) and the general method of solution is invarient. 

An additional variation would be that the substrate, rather than 
the enzyme, be elaborated from a precursor. The more general prob- 
lem is the case in which the substrate (or any Ci) ds furnished to 
the system proper by some process, physical or chemical. This prob- 
lem will be discussed in a later paper. 
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The individual cases treated in this paper include, with some ob- 
vious and clearly indicated extensions, the majority of the types of 
systems encountered in cellular respiration. The aerobic dehydro- 
genases which require either Co; or Co; would be included in the 
“typical dehydrogenase system” by merely including the steps which 
lead to the reoxidation of reduced flavoprotein. The diphosphothia- 
mine enzymes and the autoxidizable flavoprotein enzymes would be 
included in the “systems involving dissociable prosthetic groups.” In 
particular, for the Lipmann system (Lipmann, 1939) and similar sys- 
tems, one must include the association of the enzyme with Mg* and 
replace the substrate, S, by the adjunct between the substrate and 
phosphate. This latter feature may be included in the problem by 
considering the addition compound between pyruvate and phosphate 
to be in equilibrium with pyruvate and inorganic phosphate. For the 
flavoprotein enzymes which are autoxidizable, e.g., d-amino acid oxi- 
dase, it is only necessary to specify that E*- S react to give E*- H, and 
product, where E*- H, is the enzyme-reduced prosthetic group com- 
plex, and to include the reoxidation of E*-H, by O,. The steady 
state condition for [E*] is still given by (48) with k, [E*-S] re- 
placed by k, [O.] [E*-H.]. As previously mentioned the “coupled 
or branched system” can be easily made to correspond to the coenzyme- 
linked systems. The coenzyme-linked dehydrogenases may in turn be 
part of more extensive systems as, for example, in the CO, fixation 
reactions leading to the synthesis of the tricarboxylie acids (Ochoa, 
1948). It should be noted that the “coupled or branched system” 
chosen here as an example is so constituted that the parameters of 
the second, (p-consuming), system control the overall rate of the first 
system. If the second system consumed the product produced in the 
last step of the first system, the two systems would simply constitute 
a single steady state system and could be treated as such. This ex- 
ample is thus a branched system in the sense that two products 
formed in a given step suffer subsequent fates via different routes. A 
notable example of such a case is that of the metabolic reactions 
which generate energy-rich phosphate compounds and whose rates 
depend upon the phosphate acceptor systems with which they are 
coupied in a branched manner (for discussion see Lipmann, 1941). 

Finally, it should be noted that the general formulation used here 
is applicable to cyclic systems which occupy an important position in 
the present-day concept of intermediary metabolism. If the product 
B, is identical with one of the reactants, C, or C,, in the first step, 
it re-enters the system. The system is then cyclic. No new features. 
are introduced. The condition that the time rate of change of this 
component be zero is identical with the equality of the 1st and nth 
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members of the set of equations (2). This. feature, in modified form, 
has in fact already been encountered in the “typical dehydrogenase 
system” wherein the enzyme, EF, which enters in the 1st step is re- 
generated in the 2nd step. This is an excellent example of the main- 
tenance of the steady state concentration of a component which is 
conserved to the system. The requirement that [EZ] be maintained 
constant in this manner is not a mathematical overstatement of the 
problem, for this requirement is equivalent to the equality of the net 
rates of the Ist and 2nd steps. In general, this requirement will be 
equivalent to the equality of some two members of the set of equa- 
tions (2). The distinction between this physical situation and that 
giving rise to equation (43) should be carefully noted. 
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A method is given for using a set of disjoint cycles as the main body 
of a neural net that will give rise to more than one temporal response 
pattern for different afferent stimuli. The method arises out of consid- 
ering the correspondence between this type of net and certain Abelian 
groups of finite order. The consideration also gives rise to a possible 
definition of the “complexity” of this type of neural net. 


We consider here a class of neural nets consisting of a set, {Aj}, 
of afferents; a set of k& disjoint cycles, {Xi}, with m,, M2, «+: , % 
synapses respectively with the afferents; a set of internuncials con- 
nected to the cycles, but not members of them, {/;}; and lastly a set 


of efferents, {H;}, see Figure 1. 


FIGURE 1 


In the following we will use the individual cycles as our units of 
structure of any such net. We label the cycles 1, 2,----, k. We 
designate the fact that the 7 + 1st synapse in the 7th cycle receives 
an impulse at the time t + 7 + 1 by the expression 


a; (¢ a j) — (40j41 9 INj+oy °° * 5 iOn 5» 40g °°" j0t;) ° (1) 
In the right-hand side of expression (1) we have the ;a;,, as the labels 
of the synapses in the ith cycle. We will designate by t a moment 


when the only activity in the net is a set of k impulses, one in each 
of the & cycles. In general we can consider this t) as the time at 


51 


52 NEURAL NETS 


which the circuit has reached a steady state activity. For the pres- 
ent we consider only the case where one pulse only is present in any 
given cycle at any time. Also, for convenience, we will consider the 
synapse in the jth cycle which receives an impulse at t as being la- 
beled ;a,. We now define a power of a:(t) by 


aif (to) = pai (to + 7). (2) 
In addition, we define the product of two powers by 
ai (to) + as (to) = pai (to). (3) 


From expressions (1), (2), and (8) it follows that 
ai(t+q)-q=p(modn;)->-a:(t + p) =ai(t + qQ). (4) 


We shall also introduce the following notation 


IT a" (to) 
Le SS5 ops (5) 
TT Om™™ (to) 
Bi; (bo + lL) = pra; (6o'+- 1) -T; | 
i 6 
= Beinn 2 
Gi = {Bi (bo), Bi? (bo), ----, Bi™ (bo) } (7) 
Ci, ete rs i, = vf Bit (to) - Bi* (to) “+++ Bit (to) ; (8) 
Gis os Pit aCe 3-7-5 i ee cee 
(9 


L ° 
= ees hasty, 


where L = LCM (n,,----, ™%). We shall refer to expression (9) as 
the “actual product” of Gi; Gi, teeny Gi. We distinguish from this 


“actual product” another product that we call the “formal product”; 
Bis (to) e Gi + Bj!2(t.) e Gj > Gij = Gi X Gi = pt{ Bi" (to) 
» Bi (b)} . 


Theorem 1. The sets defined in equations (7) and (9) are multi- 
plicative, cyclic, growps of orders n, and LCM (ni, “sos , Ty) Tespee- 
tively. . 


(10) 


Proo f. As both are finite sets we need only show closure and com- 
mutativity under multiplication. We have 
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Bj} (to) - Bim (to) = Bii™ (to) = Bi (ty) = Bi” (to) - Bi} (to), 


from equations (2), (3), (4), and (6). The identity of the group is 
B,"(t)). This shows that the set (7) satisfies the theorem. The proof 
for the set (9) is almost identical. It will be noted that the set (7) 
is included in the set (9). The total number of groups of the form of 
set (9) obtainable from our k disjoint cycles is equal to 2*. 

The product defined in (10) can now be seen to be a direct prod- 
uct since we have all k cycles disjoint from one another. We now 
consider the general product 


Gi x Gia 2G: = Gs ests a. (11) 


Theorem 2. The formal product of 1 of the k cycles as defined in 
equation (11) is an Abelian group of order Diss 2 ty Min 

Proof. Gi; ++, 4, is the product of J cyclic groups all of which 
have the same identity ((5), (6), (7), (8), (9)). Each of the cyclic 
groups is disjoint from the others except for the identity ; hence, their 
product is direct and composes a group Gi. SPE This group is 


Abelian and its order is the product of Mi 77+ y 1. 


Corollary 1. The groups of Theorem 1 are all subgroups of the 
general group G,,».....,x- 

Corollary 2. The order of G,,2,....,, may have divisors that are 
not orders of any of the groups of Theorem 1. Suppose, for example, 
that the smallest order of the groups of Theorem 1 is 2q. Then the 
divisor 2 of the order of G,,....,, is not the order of any subgroup 
of Theorem 1. 

Corollary 3. All divisors of the order of G1,»,....,, are orders of 
subgroups of G,,....,, that are contained in Theorem 1 if, and only 
if, the following condition holds: All of the m; are primes or prime 
powers such that with any prime power all other smaller positive 
powers are represented among the 7;. * 

Proof. (a) If the n; are all primes then the only divisors of 
M, «++» %, are the primes themselves, or some product of them. All 
such products are represented as orders in the groups of Theorem 1. 
(b) If some 2», = p? and p**, pt, +++. , p are not included as the 
values of some of the other n,’s, then the ones that are missing rep- 
resent divisors that are not orders of the groups of Theorem 1 in 
all cases other than (a) and (b) of the corollary, we will have such 


divisors. 
So far we have seen that all groups included in Theorem 1 are 
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subgroups of Gy,2,....,x. The largest subgroup of G,,...,~ included in 
Theorem 1 is the group Gi,»,....,x with index equal to 


Nyse NM 
LCM (n, Dacian he) N) 
HCNF (n,,--:-, %) was defined previously (Roberts, 1948) by 


= HONF (m,---- ,%) = Ah. (12) 


k k 
HCh ( 11 m,/n: ) : 
t=1 ial 

We shall denote by a pathway the progressive excitation of sets 
of k synapses (one from each cycle) from some time t to the time 
t +p, where p is the period of the set of cycles. Thus a pathway is 
the totality of simultaneously acting sets of synapses, k-tuples, that 
can arise through a single impulse acting in all & cycles. By disjoint 
pathways (dp’s) we shall mean pathways which have no k-tuples in 
common. 


Theorem 3. If Gi,,--+* , 1, has index q in Gis “+++, 4, then there 


are q dp’s through the | cycles i,,-:--, i. 

To prove this theorem and to find the q dp’s one needs only to 
find the g — 1 cosets of Gi es ean Ga -+++,4,. Each coset and the 
subgroup itself constitute dp’s. The elements of the cosets are the 
k-tuples of the pathways. 

Characterization Matrices. In some instances the cosets of 
Gan Ge +++», 4 are not too difficult to find. We give here 


i} 
a method that depends on a matrix representation of Gis sees 


The matrix (a) of expression (13) is the matrix for Cae 
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1 1 1 i 1 1 
2 2 2 2 2 2 | 
Mi Ni eg ay 
i 0 0 0 
2 if i 1 (13) 
2 2 2 
& 
5 ‘ 0 0 0 
(b) 
Ni = N; 
1 2 q 
(a) 
If we now reduce all elements of the matrix (a) in expression (13) 
by d, where d = HCF (ni, “++, mi,), we get (b). The bracketed 
part of matrix (b) in expression (13) we call the “characterization 
matrix” (CM) of (a) or of Gi, St (If d= 1, then the CM can 


be thought of as being (a) itself. We interchange (a) and Gi, Rapes 


in our discussion as if they were identical.) If we now define the 
operation of cyclically permuting a column by one unit, that is, chang- 
ing the column {1,2,----, d —1, 0} to {2,----, d—1, 0, 1} and 
then cyclically permute an arbitrary number of columns by any num- 
ber of units, we shall get a maximum of d’ different rows (the num- 
ber of ways d different kinds of things can be put into / boxes). As 
there are d rows in the CM these d’ different rows can make up 
d'. matrices of d rows each such that all column sums are 1 + 2 + 
---+» + @d—1. The original CM is one of these d’* matrices and the 
others are corresponding CM’s for the cosets of Gi, FP Rs 
each CM it is clear that we can return to a complete matrix (a) such 
that no two CM’s give matrices (a) with an identical row.: For every 
(a) there is one, and only one, CM possible. However the converse 
is not true. =i ht 
As we have 
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Mr Ni : 
= H,, 


LOM (7; , +--+, Ni ) 


we immediately see that d'*|H, . Now we see that if d’* = H, , there 
is a 1 — 1 correspondence between the cosets of Gis ee and the 


CM’s obtained above. In this case the cosets are completely deter- 
mined by the CM’s. If, however, d’+ = H, (q an integer) we have 
a one-many correspondence between the CM’s and the cosets of 
Gi, Bernt In this case the method finds “ of the cosets uniquely , 


but the rest have to be found by trial. 

One can also merely cyclically permute the columns of (a) to get 
all of the coset representations (as one has to do when d = 1), but 
this is a more difficult task than the method of the CM’s given above, 
especially when J becomes large. Note that when d'* is equal to A, 


this means that the numbers Mis sees, ni, , obtained from Mig rrrry Ns, 


under division by d, are relatively prime in pairs. 

We have seen that two disjoint pathways are two sets of k-tuples 
of synapses such that no k-tuple is common to both sets. In some cases 
such disjoint pathways can be seen to give rise to other pathways dis- 
joint with them. As an example, consider the k-tuples (11, 12, ----, 1x) 
and (li, le, -:-+, lee, 2x1, 2.) which are members of the dp’s A., 
A, respectively and which are both acting at time t. From these two 
we see that the k-tuple (1,, 1., -:-- , 1x1, 2x) is also acting at time 
t and gives rise to another dp disjoint from A, and A.. When such 
is the case we shall say that the dp’s A, and A. are not “effectively 
disjoint pathways” (edp’s). 

In other words, any k-tuple formed in the above manner from 
the k-tuples of a set of edp’s is already a member of one of the edp’s 
of the set. 

We denote by the set of effectively disjoint pathways through a 
set of & cycles the maximum number of dp’s that can exist together 
simultaneously in a net without giving rise to other dp’s. There is 
a unique correspondence between a set of D afferents and a set of D 
efferents of the net when the net contains D edp’s. We see that the 
number of edp’s is not necessarily the same as the number of dp’s. 

If we define the quantity H. as follows: 


k k i 
H,=HCF | ( IT n;/n; +n; : (14) 
4,k=1 fmt 
i#k 
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the following three relations hold between the number, D, of edp’s 
_ and the quantities H, and H, and LCM (n,, -::: + M%)- 


HCF (n,,----,m) <D<H,, 
H,.=1>D=H,, hay gun (15) 
LOM (1, ,-++-, %) = LCM (1% ,-+:+, m1) =m? D>N. 


The exact determination of the value of D in the general case has 
not yet been made. a1 

From the above considerations we see that if we define “com- 
plexity” of a net to be the number of edp’s of the net, then the HCNF 
of the numbers of synapses in the disjoint cycles of the net will be 
an upper bound of the complexity and can be used to compare two 
nets, in a limited sense. The actual measure, ,D, of. the complexity 
under this definition is, so far, an undetermined function of the n,’s. 
If a net has a number of edp’s equal to D, then we can put D affer- 
ents into 1 — 1 correspondence with D efferents through the set of k 
cycles of the net, with each afferent giving rise to a temporal re- 
sponse pattern of any form whose period is < LCM (”,,--:-,%) =P 
and which has for the maximum number of neurons in a cross sec- 
tion taken through the cycles of the net 


k 
2k +> n; (16) 
t=1 
neurons. With the exception of a few isolated cases this method does 
require the q afferent stimuli to be in “phase” before impinging on 
the cycles. This will cause a time delay to be introduced between the 
q afferent neurons and the & cycles. 


FIGURE 2 


As an example of a simple circuit where we use the same two 
cycles for the passage of two impulses we give Figure 2 for Go, . The 
numbers in parentheses in the figure indicate the synapses in the two 
cycles respectively which send neurons to that synapse. 
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Comparing the present results to those of the previous paper 
(Roberts, 1948), we see that to each element cs ie Gis eae i 
there corresponds a set of synapses, one from each of the / cycles con- 
sidered, all acting at the time t + m. Also the “range” of the net con- 
sisting of just these J cycles is equal to the order of Gi, yin 
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